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We show how to incorporate fractionaUy charged quasielectrons in the finite quantum Hall matrix 
model. The quasielectrons emerge as combinations of BPS solitons and quasiholes in a finite matrix 
version of the noncommutative 0* theory coupled to a noncommutative Chern-Simons gauge field. 
We also discuss how to properly define the charge density in the classical matrix model, and calculate 
density profiles for droplets, quasiholes and quasielectrons. 
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CN ■ I. INTRODUCTION 



During the last few years a new class of models for the fractional quantum Hall (QH) effect has emerged. The 
I , basic construction is due to Susskind-'^, who observed that the Laughlin states at filling fraction v ~ l/(2fc + 1) are 
J~i naturally described by a noncommutative Chern-Simons (CS) theory, or equivalently, by an infinite matrix model 
c/2 \ with the lagrangian, 
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Tr{(x«-i[X",ao]m)eahX' + 20ao}, (1) 



''O ■ where X"'{t), a = 1,2, and ao{t) are hermitian matrices - the latter being a Lagrange multiplier imposing the matrix 
^ ' commutator constraint, [X-'^,X^] = i9 . The area 9 that enters Lq is the noncommutativity parameter, and B is the 
transverse magnetic field. 

As it stands, this model has only a single state, since the solution to the constraint, which can only be satisfied by 
infinite matrices, is unique (up to gauge transformations). This reflects that the theory is topological and thus has 
no excitations when defined on an infinite planei^ 

The parameter 9 can be interpreted as an area per particle, giving the unique state a constant density p — 1/2119. 
Modifying the constraint by hand, one finds other solutions correponding to fractionally charged quasielectrons and 
f — ' quasiholes^. 
\^ - In an important development, Polychronakos extended the model by supplementing |^ with the lagrangian, 

o 

^ : ifc = $t(jao_a„)$ , (2) 



where $ is a complex bosonic iV- vector— . The X^'s in QJ are now hermitian N x N matrices, and the constraint is 
changed to 



rA : [^i^X^l ==i6i-— $$^ =i6i(l-i$$1'), (3) 

O ' where k = eB9 is the so called level number. It is striking that this finite QH matrix model (QHMM) already at the 
classical level describes several key features of the quantum Hall system: 



, •In the presence of a rotationally invariant confining potential, the groundstate is a finite size circular " droplet" 

K% ' with a constant bulk density p depending on the level number. 



• The excitation spectrum is consistent with that of a QH droplet. In particular there are quasihole states in the 
bulk and gapless quasielectron - quasihole states at the edge. 

• In the absence of a potential, there is a set of degenerate low density states corresponding to single particles in 
the lowest Landau level, at well separated positions in the plane. 

In particular note that the presence of quasielectron and quasihole excitations takes this description beyond that of a 
classical incompressible fluid, and we shall see below how the model also describes how QH droplets are formed from 
well separated particles in a strong magnetic field. 

Quantizing ^, and assuming the underlying matrix degrees of freedom to be fermionic, Susskind showed that the 
density is quantized at the Laughlin fractions v = p/po = l/(2fc + 1), where k is integer and po = e_B/27r (when 



H — 1) is the density of states in a single Landau level J^ At a technical level, it was also shown in reference y that in 
the presence of a quadratic potential, there is an exact mapping of the QH matrix model onto the Calogero model, 
both in the classical and the quantum case. This mapping yields explicit expressions for both energy levels and wave 
functions^. 

In a previous paper we extended the QHMM model further, by constructing a class of conserved charges and 
accompanying currents, thus allowing for a coupling to an external electromagnetic field^. We then went on to 
calculate low momentum response functions in the classical model, in particular: 



• The ground state density, being the response to a constant electric potential, Aq. 

• The quantum Hall response an- 



• The response to a weak and slowly varying external B field. 

The results were all in agreement with the known properties of the Laughlin states. 

In spite of these successes there are several basic aspects of QH physics which are not incorporated in the finite 
matrix model given by ^ and 0. Most significantly: 

• There is no unambigous definition of density. 

• There are no quasielectron solutions. 

• There is no natural way to introduce spin and/or multilayer degrees of freedom analogous to the usual description 
in terms of multi-component CS fields^. 

• There is no generalization to fractions other than the Laughlin ones. 

In this paper we address the first two points, in some detail, while the third is addressed in another paper^. Before 
turning to the technicalities, we will give some general comments to the above list, and also briefly discuss the status 
of the quantized QH matrix model. 

In reference |y we constructed a class of conserved currents for the classical matrix model. As we will discuss below, 
there are three natural conditions on the charge density operator: it should be non-negative, satisfy the classical 
version of the sine-algebra characteristic of the lowest Landau level, and have the correct limit for particles separated 
much further than the magnetic length. Unfortunately, we have not found any definition that satisfies all these 
demands. 

The absence of quasielectrons in the noncommuative theories is related to the existence of a minimal area for each 
particle. A clue as to how to get quasielectrons is given in Susskind's original paper where they emerged from an 
ad.hoc. change of the constraint. In a recent paper, Bak et.al. showed how the addition of a noncommutatve scalar 
field (f), provides a dynamical version of this mechanism, and gives a model with soliton solutions with charge density 
larger than pQ? . In section HI we shall construct the corresponding finite matrix model. 

The problem of spin, (or pseudospin corresponding to e.g. a multilayer index) derives from the restricted nature 
of noncommuative gauge theories - U(N) is the only allowed gauge group^ifi. The standard multi-component CS 
lagrangians employed to describe spin and pseudospin, as well as the general classification of abelian QH liquids given 
by WenZ., are all based on the gauge group [/(I)'''. 

The problem of finding non-Laughlin states, as already mentioned, is superficially the same as for spin - there is 
no noncommutative version of the standard multi-component CS theories. As we show in reference's! however, the 
spin problem can be addressed by introducing fermionic degrees of freedom and couple them in a judicious way to 
the bosonic matrices. We know of no such construction for generating non-Laughlin states, and the initial hope that 
the matrix theory would provide a new and more powerful framework for the classification of QH liquids has so far 
been elusive. 

Thus, turning to quantum theory, there is no matrix model where the density is quantized to other fractions 
than the Laughlin ones (except for the trivial case of direct sums), and in particular there is no way to get the 
experimentally prominent Jain seriesi^, v — n/{2pn±l)). At a technical level the quantized QHMM is hard to handle 
since the current and density operators are mathematically very complicated objects. This means that although the 
quantum states of the model are known via the mapping from the Calogero model^, it is not possible to calculate 
density profiles. Even for the simplest case of two by two matrices the manipulation of exponentials of matrices with 
(quantum) noncommuting elements is very difficult. 

We already stressed the pros and cons of the classical matrix model, and the aim of this paper, and reference |a, 
is to extend this model to allow for quasielectrons and spin, and also to find a density operator that can describe 
quasiparticle and edge profiles consistent with what is known about the QH system. As we shall see, this endeavor has 
been rather successful, at least on a qualitative level. Within an extended finite QH matrix model, we can describe 



QH droplets, exponentially falling edges, quasihole and quasielectron excitations. In reference Ig it is also shown how 
to incorporate spin. 

The paper is organized as follows. In the next section we discuss the ambiguities in the definition of the density 
operator and give the arguments in favour of our special choice. We then calculate density profiles for droplets, and 
quasiparticles and compare with what is expected from other approaches such as CS mean field theory, and Laughlin 
wave functions. In section III we first show how to incorporate densities larger than p by adding a scalar field to 
the finite matrix model. The resulting theory has soliton solutions with integer charges, and quasielectrons can be 
constructed by adding holes on top of these solitons. We give explicit expressions for the solutions and calculate 
the density profiles which are again compared with alternative descriptions. In the last section we summarize our 
results and contrast the classical matrix model approach with the standard classical commutative CS description. 
Some technical points about the positivity of the density operator and possible alternative definitions are given in an 
appendix. 

II. PARTICLES, DROPLETS AND QUASIHOLES 

In this section we shall study the density profiles of various solutions of the finite classical matrix model. These 
solutions were all found by Polycronakos^, who also determined gross characterizations such as the radius and mean 
density of the QH droplet, and the charge of the quasihole. To calculate the profiles, we must first give a definition of 
the density operator. As stressed in the introduction, our choice, although not unique, gives profiles in good agreement 
with those obtained by other methods. 

A. The density operator 

In reference la we constructed a class of conserved currents for the classical matrix model. The general form of the 
charge and current was given by 

p{y,t) = T4{y^^ - X'^d))] (4) 

]{v,t) = Tr[(X-i[X,ao]m)^(y"-X"(t))], 

where i5(j/° — X"") is a matrix-valued kernel. The general form of this kernel follows from symmetry considerations 
and current conservation, 

6{y^ ~X^)^ j -^ fiUy'^ ~ X'^),ea,k%y'> - X% 0k') (5) 



and a rather natural guess is 



(27r)^ 



Swiy" -X")^ I ^ e^^»to"-^")g(0fc2) , (6) 



(27r) 



where g{0) ~ 1. (The special case g{x) = 1 is known as Weyl-ordering.) For almost diagonal matrices, appropriate 
for widely separated particles at positions a;„ = X„„, this corresponds to p^ = g{9k') J^n e"*'^'^". Clearly g can be 
thought of as a formfactor, and for the particular choice g{x) = e~'^^l'^^ we have gaussian "blobs" which are very 
suggestive of maximally localized one-particle wave- functions in the LLL. With this motivation we shall use 
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PI = e-^Tre-^"--^ , (7) 

where the length d was chosen as to match the rms radius \J {r^) = \/2(. = (2fi/ei3)^/^ which is the appropriate value 
for an isolated electron in the LLL. 

Before proceeding to use the formula ^ to calculate density profiles, we mention some problems with this definition. 
The first, and most severe, is that p[x) is not positive definite on the space of matrices satisfying the constraint ((JJ. 
This is not obvious, but can be shown by numerical calculations which also indicate that this is mainly a problem for 
very small systems, typically N < 10, and also gets more severe with lower v. For moderately large N, very small 
violations of positivity is seen in typical density profiles such as the "droplet" solution shown in Fig. QJfor N = 50. 
For extreme cases, such as iV = 2, the violation of positivity is large, as shown in the appendix. We have not been able 



to show that the definition ^ gives a positive definite density in the hmit N ^ oo, although our numerics appears 
to support this possibihty. 

The situation is less favourable for other ordering prescriptions. So will for instance antiordering, defined by 

Uz -Z,-z-Z^)-J-0^ e^(^-^)e*(-^^)5(^fc^) , (8) 

where again ^(0) = 1, give strongly fluctuating profiles, and large negative values for the density even for rather 
large Nm^ By going outside the class of density operators that can be written on the form Q, i.e. as a trace of 
a matrix kernel, one can define a positive definite density operator with the correct limiting behavior for separated 
particles. This construction, which essentially involves taking the square root of a delta function, has, however, other 
shortcomings. Technical details are given in the appendix. 

A second problem is that we would expect the Fourier components of the quantum mechanical density operator to 
satisfy the following commutation relation. 
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[Pk,Pp]QM = 2isin l-{k X p}\ e2 ^'Ppj:^^, (9) 

which is the sine algebra pertinent to the density operator projected onto the LLL. We have not been able to find 
any definition of the density that satisfies lO except for 9 — 0, where the anti-ordering is known to be correct. The 
claim in reference |6| that a particular quantum reordering of ^ satisfies Q for A'' = 2 is erroneousi^ Actually, by 
studying the classical limit we can show that there is no quantum ordering of neither the matrix Weyl ordered nor 
anti-ordered density operators that satisfies ©. 

For readers familiar with the string theory literature, the following comment might be of interest. In string theory 
one can show that the Weyl ordered expression for the density, corresponding to g{9k^) = 1 in ©, gives the density 
of the lower dimensional RR-charged D-branes. This follows since Q is nothing but the Seiberg-Witten map for the 
noncommutative field strength which implies that it couples to the Ramond-Ramond forms in precisely the correct 
way to act as a source of the corresponding RR-charge-'^''. In our case there is no such reason to use Weyl-ordering to 
define the density of particles and we may modify this expression as long as it respects the symmetries of the problem. 
Note however, that our choice iQ) coincides with Weyl-ordering for the k — Q component corresponding to the total 
charge. 

To summarize, we have no a priori reason to choose © rather than e.g. antiordering, or in fact any other ordering 
in the general class (0). Similarly, there is no theoretical motivation for taking any particular g(x). Instead our choice 
O is phenomenologically motivated, and its usefulness will be demonstrated in the rest of this paper. 

B. From particles to droplets 



For the gauge choice $t — \/eBO{ll ... 1), the constraint (PJ is solved by the following matrices^ 

where we (arbitrarily) chose to diagonalize the hermitian matrix X^ . For widely separated x^is, the off-diagonal 
terms, that are responsible for the " 0-repulsion" , are small, and the diagonal elements can be interpreted as the 
coordinates of the particles. More generally, we can think of the (gauge invariant) eigenvalues of the matrices Xi as 
particle coordinates Xi and j/i. Note, however, that there is no unambiguous way to pair these eigenvalues to position 
coordinates for the particles. 

Another convenient gauge choice is <I>t = \/Nk,{Q . . . 01) and introducing the dimensionless complex coordinates 
Z = -^=[X^ + iX"^), the constraint takes the form 

[Z,Z^] = l-N\N-l){N-l\, (11) 

where the bra-ket notation refers to an oscillator basis as explained in e.g. reference llOi In the large A^ limit, this is 
the usual ladder operator algebra, and the effect of the boundary field is only at the " edge" of the matrix. It is thus 
natural to seek a solution for Z similar to the lowering operator a in the n-representation. One finds, 

N-l 

Z=^V^|n-l)(n|. (12) 

n=0 



We will refer to this as the droplet solution. By a [/(]V)-transforniation it can be put on the form (|10|l . with ym — 
(since the matrix elements in Z are real) and almost equidistantly spaced x^'S, which are the (gauge invariant) 
eigenvalues of the hermitian combination {Z + Z'')/2m^ 





1 as given by the density operator ^ and the droplet solution 1121 1 with A^ 
= 2ve = 1. 



50. We use 
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FIG. 1: To the left p{x) for ly = 

dimensionless units given by 2£^ = 2v6 = 1. To the right the density p(x) in profile along the i-axis is given for // = 1, , ^ , 
respectively. All states have roughly constant bulk densities, but note the increasing wiggles, and the building up of a rim at 
the edge, for lower v. Note also the small violations in positivity at the edge. 



Using the choice ((TJ for pj^ , we can calculate the corresponding x-space density profile p{x\ which is shown in Fig. 
^for V = 1, |, -^ and N = 50. The lower v is, the more pronounced is the up-shooting rim at the edge. Excluding a 
circular segment containing the rim, the distribution is very well fitted by the formula 



.(..)^f(l-tanh^ 



(13) 



with ro « 0.99V2eN = 0.99iy/2N/v for all three i^ and /3 w 1.05, 0.92, 0.82 for z/ = 1, i, i and iV = 50 respectively. 
This is consistent with the expectation of a constant bulk density and a very rapid fall-off at the edge over a distance 
of the order of the magnetic lengthiSi 

In Fig. |21we illustrate how a droplet is formed when several, initially well separated, particles approach each other. 
The middle figure is a density plot of the droplet solution H12() for seven particles. From this solution we extracted 
the eigenvalues Xm, and then generated a set of solutions of the form (|10|) by scaling the x^'-S by a common factor, 
A. The top figure shows the density for A = 5, corresponding to particles well separated on the x-axis. In the limit 
of large A, the Xm-S are simply the coordinates of the particles. The bottom picture is for A — 1/5. Because of the 
0-repulsion the particles cannot be compressed further than the droplet, and the result is instead particles separated 
along the conjugate y-direction. That this effect is entirely due to the finite value of 9 is demonstrated in Fig. O 
which is identical to Fig. |21 but with the off-diagonal 0-repulsion terms in IjlOII set to zero (still keeping the same 
formfactor g{9k'^). In this case no circular droplet is formed and the maximally compressed state is simply an overlap 
of the individual gaussian distributions. 



C. The quasihole solution 



The droplet solution (|12|l can readily be modified to describe a quasihole, i.e. a state where the density close to 
the origin is depleted compared to the droplet state. Polychronakos found 



N-l 



Z = y?|iV)(0| + Y, V^^+q\n - l)(n| , 



(14) 



where < q corresponding to a shift in the eigenvalues of the radius operator 

N-l 

R' = {X'f + {X^f = 2eY,{n+-+ q)\n){n\ , 



(15) 



n=0 



with the amount q relative to the original droplet. By inserting H14|) into Q, we get the distribution p{x) shown in 
Fig. H 
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FIG. 2: The compression of seven quantum Hall particles as discusssed in the text. Note the perfect circular symmetry of the 
maximally compressed state in the middle picture. 
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FIG. 3: As in FigElbut with 6 = 0. No droplet is formed and the maximally compressed state in the last picture is just seven 
superimposed gaussians. 




FIG. 4: The density distribution p{x) and its profile given by Q for H14^ for A^ = 50, u = q = 1/3. 





FIG. 5: To the left, Monte Carlo calculation^'^ of the difference Qhoic(-R) — Qdropict(^) in the cumulative integrated charge 
between the Laughlin quasihole density distribution and the ground state density distribution for A^ = 50, v — 1/3. The 
horizontal line corresponds to the charge difference 1/3. The radial distance is in units of \/^£ and the charge in units of e. To 
the right, the same distribution calculated from the matrix model. 



Although there is a clear charge deficit at the origin, the matrix model does not reproduce the complete expulsion 
of the electrons characteristic of the Laughlin quasiholes. A more detailed comparison is made in Fig. |5l where 
we show the cumulative integrated charge Q{R) = tt L dr^ p{r) for the Laughlin quasihole (left)i& and the matrix 
model quasihole (|14|l (right) . We also calculated the root mean square radius for a quasihole of charge ve in a state 
with filling fraction v numerically with the results \f(r^ jl — 2.6, 1.4, 1.1 for v = li |j ■! respectively. This can be 
compared with the vortex solution of the mean field composite boson model (see e.g. reference il7i) . where the vortex 



has 



= v2.he. 



III. BPS SOLITONS IN THE FINITE MATRIX MODEL 

We already mentioned that the finite matrix model defined by ^ and Q docs not allow for quasielectron solutions. 
On the other hand, such solutions can be found if the constraint is modified by hand. In the infinite matrix model 
we can take 



[Z,Zt] = l + q|0)(0|, 



(16) 



which describes a quasihole at the origin for g > and a quasielectron for q < 0. To have dynamical quasielectrons a 
constraint of this type has to appear as one of the equations of motion. Such a construction, based on a noncommuative 
version of the Chern-Simons-Higgs model, was given by Bak et.al. in reference |9|. We first briefiy review their work, 
and then show how to construct a corresponding finite matrix model. This will require both a modification of the 
action for the noncommutative scalar field, i/i, and a coupling between <^ and the boundary field $. 

In terms of the complex covariant position operators Z = — i=(A'^ + iAT^) and Z^, and the CS level number k, the 

noncommutative CS lagrangian J^l takes the form. 



Lcs - yTr {Z^^^Z - ZDqZ^) + nTiao, 



(17) 



to which Bak et.al. added a noncommutative 0^ lagrangian. 
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(18) 



Here is a matrix field in the fundamental representation, i.e. it transforms as ^ -^- U(j) under the gauge transformation 
U . The covariant derivatives are defined by D^^ — d^ 



ia^ and act on as 



DqcJ) = docf) + ia,o(f> 

1 



-e,j[Xj0-0% 



(19) 



where Xi are noncommuting coordinates, [a;i,a;j] = itijO. The corresponding derivatives are given by di = ■z^ij[ 
and the matrices Xi, defining the actual state, are related to the noncommutative gaugepotential via 

VVo — Xj U^inCln , 



(20) 



i.e. hi parametrizes the deviation from the ground state solution [X^, AT^] = iO. 



Defining the current 



■h = ^[(D,<^)0t_^(D,<^)t], (21) 



ttie Hamiltonian can after some algebra be written as 

H = Tr-^(I?i0±zi?2</')(i?i</'±ii?2</')^±^TrAJj±7r^Tr([Z,Zt]-lTAm#t)^^t. (22) 

Assuming the solution to be regular enough for the covariant derivative of the current to integrate to zero, and taking 
Xmn = ±1 so that the last term will vanish because of the constraint 

[Z,Z^] = l--c^<j,\ (23) 

K 

the Hamiltonian reduces to the first term which is quadratic and equals zero for {Di(j) ± iD2<j)) = 0. Thus this choice 
of parameters corresponds to the theory being of the Bogomol'nyi-Prasad-Sommcrfield (BPS) formSiiS. The complete 
set of BPS equations is given by, 

ao = ^^H^ (24) 

Dic/)±iD2<l) = 

[Z,Z^ - l--00t, 

K 

and it is easy to check that any solution of these is also a solution to the full time-independent equations of motion 
corresponding to Lcs + ^s- 

We now turn to finite matrices. Because of the coupling ^ Ti{aQ(j)(j)^) the Gauss law constraint became (|23ll . which 
can be satisfied by finite matrices. Of course, it is then no longer possible to have [x^,x^] = i6[z,z''] — i9, where 
z = -^{x^ + ix"^)- Instead we let 

[z,z'']= M =1- N\N-1){N-1\, (25) 

with 

JV-l 

5=5^V^|n-l)(n|. (26) 

n=l 

Now the Hamiltonian can no longer be written on BPS form l|22() . but by adding the term, 

Lm - T^Tr(+0A/0t _ 0^t) ^ (27) 

it is a matter of algebraic manipulations to show that the Hamiltonian corresponding to L = Lcs + Lg + Lm is again 
of the BPS type. 

It is not hard to verify that the model we just defined has droplet solutions, and topological solitons of the type 
found by Bak et.al?. There are however no quasihole solutions. This can be remedied by also adding a Polychronakos 
type boundary field, which has the additional advantage that the sector where the scalar field is not excited becomes 
identical to the original finite QH matrix model. Our final lagrangian now reads 

L = Lcs + Ls + Lm + L^ (28) 

where the boundary lagrangian is given by, 

L$ = $^iDo$ - ^$^0^* (29) 

2c' 

yielding the Gauss law constraint 

[Z, Z^] = 1 - i(/)0t - l$$t _ (30) 



The last term in H29() was added to allow the Hamiltoninan to have a BPS form almost identical to (|22() . but with the 
BPS equation H23|) replaced by (|30|l . The remaining BPS equations are unchanged. This completes the derivation of 
the extended finite QH matrix model, which is a finite matrix version of the conformal Chern-Simons-Higgs model 
introduced by Jackiw and P:^^. 

We now turn to a discussion of the solutions of this model. First note that all solutions discussed in section II, 
i.e. the isolated particles (|l()|l . the droplet (|12|l . and the quasihole (|14(l . can all be taken over unchanged if we set 
= 0. For non-zero (j) we will have two new types of solutions corresponding to waves and solitons. The latter, which 
will provide the basic building block for the quasielectrons, are the most interesting, but we first briefly discuss the 
former. 

A. Collective modes 

For our model to give a realistic description of the QH system it is important that the collective wave-like solutions 
in the bulk are gapped. This is certainly expected from the analogy with the continuum model, but should nevertheless 
be established in the matrix model context. Let us first consider the case of a constant density of (p particles, p, (not 
to be confused with the constant density of electrons p = -^ represented by the solution Z — z,(j) ^ 0) in the infinite 
matrix model of Bak et.al.^. The mean field solution that we want to expand about is given as an expansion in the 
density p 



rp{i + o{p)) (31) 



ao = ~^p + Oip^) 

which solves the full equations of motion to first order in p. 

We then expand ((T7|l and ((TH|) to quadratic order around the mean field solution, and then use the polar decom- 
position (j) = UP which is valid for an arbitrary square matrix. Here U is a unitary, and P is a positive semi-definite 
hermitian matrix-'^^. By a U{N) gauge transformation, we can now remove the U{N) phase U from the field (j). As 
a result we find that the kinetic term of what remains of becomes a total derivative and this field thus becomes a 
Lagrange multiplier enforcing a constraint relating the fluctuations of the gauge field and cj). The result is that we 
have moved the entire dynamics from the scalar field to the gauge field Z. This is the noncommutative version of 
going to unitary gauge in Ginzburg-Landau Chern-Simons (GLCS) theory. The resulting lagrangian reads, 

Lcs + La, = inTia'a — - — Tr(aa^ + a^a) + . . . , (32) 

in terms of a = i{z— Z) = ^J9 /2{ai +ia2)^ with ai given by H2Q(I . The dots indicate commutator terms corresponding 
to spatial derivatives, as well as potential terms and terms of higher order in p. There is also a constraint equation 
that relate density fluctuations to the noncommutative gauge field. Just as in the commutative case, the lagrangian 
(|32|1 has the form of a harmonic oscillator, and consequently exhibits a gap at lOc — pj nm. This has the natural 
interpretation as the Kohn mode at the cyclotron frequency of the </> particles. 

For solutions with vanishing background density p - the simplest case being that of a soliton considered below - 
there will be zero modes corresponding to translations. In the full model (|28|) there will also be gapless edge modes. 
We have not analyzed these more complicated cases, but we think that the above demonstration of the similarity 
between the commutative and noncommutative models strongly suggests that the latter will not develop any gapless 
modes not found in the former. 

B. Solitons and quasielectrons 

In reference |9|, Bak et. al. found noncommutative counterparts of the self dual vortex solutions due to Jackiw and 
Pi. These correspond to a quantized flux, and, as will be clear from the explicit expressions given below, they carry 
unit electric charge. Since flux is quantized, one cannot have fractionally charged quasielectrons in the model by Bak 
et.ai, but in our finite matrix model there is a natural construction in terms of a soliton combined with a quasihole. 
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FIG. 6: The density distribution of a soliton of charge —1 for A'^ — 50, v = 1. 



Using p9(l one can derive 



1. The charge —1 soliton 



Di4)- iD2(i 



;(^v-#^): 



(33) 



and we will take the sign in H22|) such that this is one of the BPS equations (|24|l . A soliton of charge —1 centered at 
the origin is now given by the following expression, 



N-l 



Z = Yl V^n^^\n- !)(?■ 



$ = \/iN - l)K|iV- 1) 
= V^|0)(0|. 



(34) 



It is easy to verify by direct substitution that this indeed is a solution to the BPS equations corresponding to the full 
QHMM 1211) ■ In Fig. Owe show a density plot of this solution for N = 50. 



2. The quasielectron 



In analogy with the soliton solution (|34|l . we can now try to construct a fractional quasi-clcctron based on the 
Ansatz, Z = J2n=i V" + 'Zl'^ ^ l)("^l with — 1 < g < 0. Note, however, that this implies Z\Q) 7^ so the only 
option seems to be = ^/(W~\~'q)K\N — 1)(0| and $ = y^—qK\0). Such a solution would however not reduce to the 
ground-state H12I) as q ^ 0. 

An obvious alternative construction, alluded to above, is to add a quasihole to the soliton, i.e. to combine the 
solutions H34I) and (|14ll . This amounts to finding a g-dependent modification of Z for which the constraint is still 
q- independent. The solution is 



iV-l 



Z - ^\N-l){l\+J2V^^^\n-l){n\ 



n=l 



$ = Vn^\n - 1) 
^ - V^|o)(o|, 



(35) 
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FIG. 7: The density distribution of a charge — ^ quasielectron composed of a charge +| quasihole on top of a charge -1 soHton 
for iV = 50, f = 1/3. 





FIG. 8: The same plots as in Fig. |3 but for a Jain quasielectron«lS> Note the characteristic dip at the foot to the right of the 
quasielectron. 

for which we still have Z\0) = and hence I)(j> = 0. Again the BPS equations can be verified by direct substitution. 
In FigHwe show for N — 50, i' — 1/3 a. q — 2/3 hole sitting on the top of a soliton of charge —1 to produce a 
charge —1/3 quasielectron. The small dip at the foot of the peak has also been seen in numerical studies of QH wave 
functionsi^, see Fig. |H1 

Note that this solution does not reduce to the ground state as g ^ 1. This is not necessarily a drawback, since 
it might be interpreted as the limiting case of a small exciton, i.e. an overlapping state of an electron and a hole. 
Finally we should mention that we have not investigated the stability of our quasielectron solution, so we cannot be 
sure that it will not decay into a soliton and a quasihole. Although such a calculation amounts to a straightforward 
small oscillation analysis, it is algebraically complicated, and also of limited interest since we would in any case have 
a stabilizing Coulomb interaction in a more detailed model. 



IV. SUMMARY AND DISCUSSION 



To summarize, we have argued for a particular expression for the charge density in the classical QH matrix models, 
and shown that with this definition, the various solutions corresponding to separated particles, droplets and quasiholes 
are reproduced in reasonable agreement with standard treatments based on wave functions and GLCS mean field 
solutions. We furthermore extended the model to incorporate densities higher than that of the groundstate, and found 
quasielectron solutions. Again the profiles were in good agreement with those found from explicit wave functions. 

In this connection it is fair to ask what has been gained by the classical QH matrix model as compared to the usual 
mean field CS theories, so we shall now briefly contrast these approaches. That the two theories are closely connected 
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is clear from Susskind's original formulation of the QH matrix model as a noncommutative CS theory described by, 

£ = - — e'""^ ( a^ * d^ax + —a^ • a^ • oa j , (36) 

where the Moyal star product • is defined with a noncommutative parameter 9 — (27rp)~^. This is a purely topological 
theory, consistent with the infinite matrix model having a unique state with a constant density p. The corresponding 
commutative CS theory is also topological and has a unique state on an infinite plane. Here we should note that this 
commutative effective CS theory can be derived from the GLCS theory by expanding around a mean field:^-^ 

Quasielectrons and quasiholes can be introduced by hand in the infinite matrix model by changing the constraint. 
In the CS theory this correspoinds to adding delta function sources. Here we see the first advantage of the matrix 
model in that it gives a size ~ 6* to the quasi particles. ^^ 

Adding the boundary field $ to the matrix model allows for a plethora of states not described by the usual CS 
approach. Defining the latter on a manifold with a boundary gives edge degrees of freedom corresponding to chiral 
Luttinger liquids, but there are no excitations inside the bulk nor outside the droplet. The basic reason is that the 
edge excitations in the CS theory can be understood as hydrodynamic modes of an incompressible liquid, while the 
matrix model allows for density fluctuations in the fluid itself. From this it is also clear that no questions regarding 
density profiles or effective sizes of quasiparticles can be addressed in the framework of pure CS theory. 

There is an asymmetry between quasielectrons and quasiholes in the matrix model, since there is a maximal density 
given by the noncommutative parameter 9. This was the basic reason that forced us to introduce a new field, (/>, 
to describe quasielectrons, while quasiholes were present already in the model based on only a CS field X^ and the 
boundary field $ needed to "absorb" the anomaly. Such an asymmetry is present also in other descriptions of the 
QH effect. For instance, in the wave function approach Laughlin's quasihole wave function is essentially unique, 
while there are several quite different approaches to the quasielectron state^^. The introduction of a new field raises 
questions about the correct counting of degrees of freedom. The finite matrix model without any extra field describes 
N particles, but with a phase space repulsion giving a maximum density ^ 1/9. As we have shown, the extra field 
relaxes the maximum density constraint in a way consistent with QH phenomenology, but one might worry that we 
have at the same time introduced additional unphysical (gapped) excitations in the high energy part of the spectrum. 
We have not investigated this problem any further. 

To summarize, there are some aspects of QH physics that is more easily described in conventional CS framework, 
notably the classification of abelian QH liquids developed by Wen^. The QH matrix model, on the other hand, allows 
for a more detailed analysis of density profiles and a dynamical description of quasielectrons and quasiholes. 

There are several detailed questions left open concerning the details of the classical QHMM, and the quantum 
theory is to a great extent unexplored territory. If we might venture a guess, we would however say that if the 
noncommutative approach to QH physics is to provide any essential new physical insights one has either to find ways 
to generalize the quantum models - with the aim of understanding the hierarchy and/or the Jain states - or to find 
some quantitative use for the classical description. The mere fact that a classical model can do so well in describing 
a strongly interacting system in the extreme quantum regime is in itself intriguing, and it might be quite interesting 
to extend the model to include disorder and study possible phase transitions. 

Acknowledgment: We thank Alexis Polychronakos for interesting discussions and helpful comments on the 
manuscript. THH and AK were supported by the Swedish Research Council and the work of RvU was supported by 
the Czech ministry of education under contract no. 143100006. 

APPENDIX A: MORE ON DENSITY OPERATORS 

In this appendix we first demonstrate that the definition of p{x) implied by |(7J) gives negative values for certain 
configurations satisfying the constraint l|2Il. We then give two alternative definitions of p(x) which are positive, but 
have other difficulties. 

1. The Weyl-ordered density is not positive 

For small electron numbers and small filling fractions one can find many solutions for where the density becomes 
negative. In some cases it just about becomes negative, but for other solutions the violation of positivity is big, as 
can be seen in Fig. |^ 

That the density is sometimes negative can, for special cases, also be established analytically to lowest order in 0. 
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FIG. 9: Density profiles of tiie droplet solution 11211 for N = 2, p = 1/5. The definition ||7|l is given by the broken line and the 
positive definite definition 1A7II by the full line. Both definitions give a circularly symmetric distribution. The density 1A7II 
integrates to a total particle number of about 1.66. 



2. Alternative density operators 



We now give two alternative definitions of the density operator that are both non-negative. The starting point is 
the non-relativistic density operator for N point particles, 



N 



p{x) = '^6{x-Xn) ■ 



(Al) 



If Xn are taken as quantum operators, this is also the first quantized density operator in the x representation. The 
operator ljAl|l is by construction non-negative, since it is a sum of positive operators. In momentum space, the 
operator IjAip takes the form 



JV 



p(^) = E' 



-ik-Xn 



(A2) 



It is well known that quantum mechanical particles in the lowest Landau level are described by the following density 
operator. 



N 



PL 



(fc) = Ee-*^""e-*' 



(A3) 



with Zn = xl^ + ix"^ and [zm, ^n] = 2€i5mn- In X space this becomes. 



Pl{x) 



1 ^ 

n=l 



N 



= '^5fXx~Xn) 



(A4) 



where 5(, can be thought of as a regularized delta function. Again the operator l)A4|l is positive by construction. 

With these preliminaries, we now present two possible definitions of p{x) in the finite matrix model that are 
manifestly positive. The most obvious idea is to try to extract a set of N particle positions, Xn from the matrices X' 
and simply plug these into a formula of the type (|A4|I . In this case we are of course free to use any positive definite 
profile function for the particles, but by choosing exactly ljA4|l we ensure that the profile of a single particle in the 
matrix model is identical to that of an electron in the lowest Landau level. The problem of defining coordinates in 
the QH matrix model was discussed in a paper by Karabali and Sakita^^. They showed that taking the eigenvalues 
of the complex matrix Z as particle positions;^ correctly reproduced the low momentum part of the Laughlin wave 
function, while the short distance part was distorted - the characteristic \zi — Zj\'^'^ behaviour of the two particel 



correlation was softened to a lower power. We would thus expect that a density operator defined by (|A4I) and the 
coordinates proposed in reference l22 in spite of being positive, would have difficulties in describing the profiles studied 
in this paper, which vary rapidly on the order of a magnetic length. Since the construction is very indirect, we also 
do not have any closed expression for the density and current similar to 1^)1. 

Another possibility is based on expressing the operator HA4|) as a square of an operator, thus making the positivity 
manifest: 



Pl{x) = {pL{x)^f 



(A5) 
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where 



N , N 



Pi(af)^«-— =^e-^i^«^(5|(f-f„), (A6) 

n— 1 n— 1 

which is a good approximation when the particles are far apart. Going to Fourier space, where the square of the 
distribution become a convohition integral, we are led to the following proposal for the density operator. 



2^2 



pp^^(p) ^ f^ / d^ke''^^P-'^r+k']T,[e-^iP-k)■x]rJ.^^^-^k■x^ (A7) 



The corresponding ppos{x) is positive by construction, and it is easy to show that for widely separated particles, where 
the matrices become almost diagonal, the profile reproduces the one given by ljA5|l and ljA6p . When the particles 
come closer this is no longer true. Figure |51 shows the droplet solution with the old definition (jT)) shown by a broken 
line, and the definition HA7|I by a solid line. Other examples, like for particles further apart, show again that if the 
density becomes negative, ljA7|l repairs that. The basic problem of that definition is that it is not normalized, i.e. 
J (Px Ppos{x) ^ A^.This can of course be remedied by a renormalization, but difficulties remain. 

Note that the definition (|A7|) is not in the general class l@J since it involves the product of two traces rather than a 
single trace over a matrix kernel. This in particular means that our construction of a conserved current is no longer 
valid, but more importantly, that Pandora's box is opened - why should we restrict ourselves to the product of two 
traces? Why not several, or perhaps even an infinite series? 

In summary, we have given alternative constructions of the density operator which are manifestly non-negative. 
There are however other difficulties related to these proposals, and we have no reason to believe that they would 
provide a better description than ^ that we used in the main text of the paper. 
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in the quantum matrix model the relation between density and level number is shifted toi/=l/(l + K) because of quantum 

fluctuations. There is however also a correction to the relation between the level number and statistics that was used in 

reference Q. These effects in fact cancel and leave the original result of Susskind unchanged. A topological argument for the 

quantization of k to any integer was given by Nair and Polychronakos-. 

The error is corrected in arXiv:cond-mat/0304271 /2. 

These eigenvalues are actually the roots of the A''-th order Hermite polynomial as pointed out by Polychronakos in reference 

u 

Taking the logarithm of the density shows that the tail in fact falls faster than the exponential in equation 1131 . A better 

fit is given by e~^^'^~'^°''^ ' . However, fitting of just the tail is much more uncertain than is a fitting where also the plateau 

is used . 

Since the quantum GLCS theory is a direct rewriting of the microscopic theory (using a singular gauge transformation) , it 

is an interesting open question whether the noncommutative action^" could be obtained by a more sophisticated mean field 

approach. 

By introducing further terms in the expansion around the mean field in the commutative CS theory, one generates terms 

~ 6^, where b is the CS magnetic field. Such a term will give a size to the quasiparticle, as discussed in e.g. reference l2ll 

All but a set of measure zero of the complex matrices Z can be diagonalized as Z — XEX~^ . 



